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It has been proposed that four-dimensional QCD with fermions in the adjoint representation ex- 
hibits volume-independence in the laige-N hmit. If correct, this would mean that results for phys- 
ical quantities could be obtained from the single-site version of the model. A necessary condition 
for volume-independence is that the {Zn)"^ center-symmetry of the single-site theory is unbroken. 
We explore the phase diagram of the theory with a single Dirac fermion using Wilson fermions 
for a number of colors in the range 8 < N < 15, and identify the region in the parameter space 
of quark mass and gauge coupling where the symmetry appears to be unbroken. Our evidence 
suggests that this region includes both light and heavy quarks, and our results are consistent with 
this region extending to the continuum limit. 
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1. Introduction 

Long ago, Eguchi and Kawai (EK) proposed the striking idea that the properties of gauge 
theories become independent of the volume when the number of colors, A^, is sent to infinity 
Using a lattice regularization, this would allow one to determine physical, infinite-volume, quan- 
tities from calculations on a single lattice site, with both color and space-time degrees of freedom 
packed into the four link matrices. This idea is of both theoretical interest — can such matrix models 
indeed represent field theories? — and of possible practical utility — it might allow the determina- 
tion of laige-N properties in a more computationally efficient way. This would in turn allow one to 
connect to approximate analytic approaches to laige-N QCD, as well as to learn about QCD itself. 

The application of this "volume-reduction" to pure gauge theories turned out not to be sim- 
ple. As discussed below, the center symmetry of the SU {N) volume-reduced lattice theory 
must be unbroken for the equivalence to infinite volume to hold, but this symmetry is, in fact, 
spontaneously broken by the single-site theory [Q, |3|]. Various proposed "fixes" — in particular, 
"momentum-quenching" and twisting — have turned out not to work. See [Qj for references to this 
history. What does work, however, is to reduce to a lattice whose physical extent exceeds a scale 
of about 1 fm [||], and this approach, which, for a lattice spacing a ^ 0.1 5fm, allows one in practice 
to work on lattices of size ~ 6^, has been successful in determining several large- volume physical 
quantities. 

In our present work (summarized here but described in much more detail in Ref. [^) we return 
to the possibility of complete volume reduction to a single-site theory, but in the context of QCD- 
like theories with fermions in two-index representations. In particular, our evidence suggests that 
such reduction holds for the theory with "quarks" in the adjoint representation. In addition, it 
appears that these quarks can be heavy, with m S> Aqcd> in which case the long distance theory is 
the pure-gauge theory. In that case the reduction we find implements the original idea of EK. 

Our motivation comes from the following relations between theories. Start with QCD, with 
N = 3, but consider the fundamental representation Dirac fermions i2Nf of them) to be in the 
equivalent antisymmetric (AS) representation with two antiquark indices. Now take N ^ °° and 
obtain the "orientifold" \aige-N limit of QCD — one in which quark-loops are not suppressed. 
This theory is "within 1 /N" of physical QCD. Next use the "orientifold" equivalence of Armoni, 
Shifman and Veneziano [j8|], which states that the even C sector of the laige,-N theory with 2Nf 
AS Dirac fermions is equivalent to the corresponding C-even sector of the large-N theory with Nf 
Dirac fermions in the adjoint irrep. Finally, follow Kovtun, Unsal and Yaffe and use an orbifold 
equivalence to relate, when N oo^ the latter field theory to the single-site gauge theory with Nf 
adjoint Dirac fermions [^. If these equivalences hold, the large-N single-site theory gives results 
for two interesting field theories, one of which is close to physical QCD.^ 

In our work to date we have taken Nf = 1 — so that the corresponding QCD-like theory has two 
degenerate quarks — but Nf = 1/2 and 2 are also of considerable interest (the former is related to 
SUSY Yang-Mills and to 1 -flavor QCD, the latter to the potentially nearly-conformal theory with 
2 adjoint Dirac fermions). 

' For details on how one extracts long-distance quantities such as the string tension and meson masses from a single- 
site model see Refs. |[l[|[ ^ ||. 
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A crucial question is, of course, what conditions must be satisfied for the equivalences to hold. 



The necessary and sufficient conditions are as follows | |10| ] : 

1. Large-N factorization holds. 

2. The orientifold equivalence requires that C is not spontaneously broken in either the theory 
with AS or that with Adjoint fermions, at large volumes. 

3. The orbifold equivalence requires that translation invariance is unbroken in the large-volume 
theory with Adjoint fermions. 

4. The orbifold equivalence also requires that the center symmetry of the single-site theory 
with Adjoint fermions be unbroken. 

The first three conditions are expected to hold (although C can be broken for small volumes, and 
translation invariance can break at nonzero density) and it is the last condition that is the weakest 
link. Indeed, as noted above, this is the condition that fails in the pure gauge case. Thus our initial 
focus has been on determining whether the last condition holds. 

There are two previous calculations which suggest that the center symmetry of the single-site 
theory may be unbroken. Ref. calculated the one-loop effective potential for Polyakov loops 
in continuum regularization on x S^, i.e. with one direction compactified. A weak-coupling 
calculation is then justified if the compact direction is small enough. It was found that, if one uses 
periodic boundary conditions, and for massless fermions, the Z/v center symmetry is unbroken for 



Nf = \ /2, 1 and 2, but broken for Nf = 0. (For a recent extension to nonzero masses see Ref. [11].) 
In words, the fermionic contributions to Veff lead to a repulsion between the eigenvalues of the 
Polyakov loop that overcome the attraction caused by the gluonic contributions. The same has been 
found to be true in a lattice regularized calculation [p^], although there remains some uncertainty 
over what happens if only the compact direction is discretized and when nontrivial forms of center 



symmetry breakings are studied 013|]. These results are encouraging, but do not directly apply 
to the single-site model. Perturbation theory for this theory has severe IR divergences, and the 
appropriate analysis has not been done. For this reason, and because we are interested in practice 
in intermediate values of the coupling where one-loop perturbation theory is not accurate (b ~ 0.33, 
corresponding to j8 ~ 6 for A/^ = 3), we have chosen to do a non-perturbative calculation. 

2. The (possibly) equivalent theories 

We choose to use Wilson fermions and the Wilson gauge action. The (large-volume) lattice 
field theory then has the action 

Sft =2NbY, ReTri7piaq + yDw V (2-1) 

plaq 

where Upiaq is composed, as usual, of the fundamental irrep gauge matrices, b = {g^N)^^ is the 
inverse 't Hooft coupling, and the fermion operator is (x and y labeling lattice sites) 



1(1- 7,) uS 5y,+, + (1 + 7,) US^ Sy,-, 



(2.2) 



3 



Volume independence oflarge-N QCD with adjoint fermions 



Stephen R. Sharpe 



with the Unks here being in the adjoint representation. There is a single adjoint Dirac fermion field, 
corresponding to Nf = I. Boundary conditions are periodic in all four directions on all fields. 

Notable symmetries of this theory are the usual center symmetry and an S0{2) flavor 
symmetry. The latter can be seen by writing the action in terms of two Majorana fields. 

The corresponding single-site theory is obtained simply by removing the sum over sites. One 
then has a theory of four link matrices and a single xj/ and i//. The fermion operator becomes 

4 

(2.3) 



This looks peculiar because the first derivative part connects back (necessarily) to the same site, 
and so vanishes ifU^ = I. This part becomes non-trivial, however, since in the reduced model the 
Un fluctuate around nontrivial background gauge fields, and one can see in perturbation theory how 
space-time (or, more precisely, momentum space) becomes embedded in the link matrices. 

The symmetries of the single-site theory are an odd-looking gauge symmetry [U^ —^D.U^Q.^ 
for all /I with Q. G SU (N)] and the main player in the following — the center symmetry: 

f/Ai^^AiZ"" with z = e'^'"^^ and n^eZN- (2.4) 

It is important that the putative equivalence applies when both theories having the same bare 
couplings b and K. The continuum limit is to be taken after N —^ooso that (if the equivalence holds) 
the single-site theory has given one information about the infinite volume field theory. 

In order to test the equivalence, it is useful to have some knowledge of the phase diagram in 
the K — b plane of the large-volume \aige-N adjoint fermion theory. Since, to our knowledge, this 
theory has not been simulated previously for any N, we begin by giving an educated conjecture for 
the phase diagram. This is shown in the left panel of Fig. |l|. The continuum limit is at the top while 
the pure gauge theory lies along the left axis. 

The most important feature is the line, which represents the line (or region of Aoki phase) 
along which (within which) the quarks attain their minimum mass. The presence of this line/region 
is predicted by chiral perturbation theory including lattice artefacts, for which the analysis is very 
similar to that for QCD, except that the chiral symmetry breaking pattern is here SU (2) S0{2). 
For more details see Ref. [^. Just as in QCD, one finds that, near the continuum Umit, one can 
either have a first-order or Aoki-phase scenario. Our data (discussed below) indicate that it is 
the former which occurs. For very small b, however, the strong-coupling expansion predicts an 
Aoki-phase. Thus we have drawn the line becoming a region at small b. For large b the line must 
approach = 1/8. 

The other feature is the (possible) line of bulk transitions running roughly horizontally. We 
know from simulations that this transition occurs in the \wge.-N pure gauge theory, and that it is 
strongly first order. It is not a deconfinement transition, but rather a lattice artefact. This transition 
must extend for some distance into the plane, but we do not know how far. Our numerical evidence 
turns out to be inconclusive as to whether the transition turns into a crossover as one moves into 
the plane. In either case, we want to work on the continuum (large b) side of the transition region 
to avoid lattice artefacts. 

The right panel of Fig. |l| shows our conjecture for the phase diagram of the single-site model. 
In this case we must take N ^ °° for there to be true transitions, although it turns out that this is 
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Figure 1: Conjectured large-A^ phase diagrams for theories with a single adjoint Dirac fermion. Left panel: 
infinite volume; right panel: single-site. 



not much of an issue in practice. The left axis (at k" = 0) corresponds to the EK model, which has 
the above-noted center symmetry breaking transition (really, a sequence of transitions) at ~ 0.19. 
Based on the perturbative calculations of Ref. [12] we expect the symmetry to be restored when 
the quark becomes light, and that it remains unbroken for a range of k: including jq.. Our con- 
jecture therefore contains a funnel-shaped region of unbroken center symmetry in which volume- 
independence should hold. We recall that, in the continuum limit, theories equidistant from on 
the two sides are physically equivalent, so we only need to work on one side. 



3. Some numerical results 



We have simulated the single-site model using the Metropolis algorithm and calculating the 
determinant (which is real and positive) exactly. CPU time scales as A^*^ for a complete update of 
each link matrix, and we have been able to work only up to A/^ = 15 using desktop PCs. Never- 
theless, this appears to be sufficient to map out the phase diagram, since the results do not change 
substantially between N = S and N = 15. For more details, see Ref. [^. 

We have tested our conjecture for the phase diagram by doing horizontal hysteresis scans for 
couplings in the range = 0.1 — 1, vertical scans for a variety of values of K, and a few high 
statistics runs at selected points in the putative funnel region. Very long correlation times limit 
our simulations to b < 1. We note, however, that = 1 is a very weak coupling, corresponding to 
/3 = 18in SU{3). 

We calculate the average plaquette as well as a number of order-parameter for symmetry- 
breaking: Polyakov loops = Tr{U^)/N and "corner loops" M^v = ^^{U^Uv)/N. In the momentum- 
quenched EK model we found that one needed the M^v to uncover the symmetry -breaking [Q], and 
we find the same here: there are regions of symmetry breaking in the K> Kc part of the phase- 
plane for which the P^i vanish but the M^y do not. To make sure that we do not miss any signals 
of symmetry breaking, we have, for a few points in the plane, also calculated the 14641 different 
traces 

Kji = Tr{U'l'U"W"'U2') /N , with = 0,±1, ±2, . . . ,±5 , (3.1) 
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Figure 2: Scatter plots in the complex plane for all four Polyakov loops for = 10, b = 0.5 and (moving 
from left to right) K" = 0.0001, 0.03, 0.06 and 0.16. Both horizontal and vertical scales run from — 1 to 1. 
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Figure 3: Scans of the plaquette as a function of if for = 8 at = 0.1, 0.18, 0.26, 0.34, 0.42, 0.50 and 1. 



where <0 implies hermitian conjugation. 

Our results are consistent with the conjecture described above; see Ref. |j^ for complete de- 
tails. As an example, we show in Fig. § scatter plots of Polyakov loops for N = 10 and b = 0.5. For 
K" = 0, symmetry breaking is clear: all four end up proportional to g2;tm^/io (-^jfj^^ jj^ {j^j^ case, 
different n^), although, for this small A'^, there is still some tunneling between different "vacua". 
Vestiges of symmetry breaking are still apparent at K" = 0.03, but are absent at K" = 0.06. For this 
value and for K" = 0.16 (just below K^) the fluctuate around the origin and average to zero within 
errors. Results atN = 15 are consistent, with the expected reduction in tunneUng for small K (none 
is seen) and the expected reduction in the size of fluctuations. 

Evidence for the presence of a first-order transition is shown in Fig. ^: there is a jump in the 
average plaquette for each value of b. We also observe hysteresis at this "transition", and a con- 
sistent behavior for A/^ = 8 — 15. The transition is not, however, associated with center symmetry 
breaking and we interpret it as an example of the first-order scenario predicted by chiral perturba- 
tion theory. This is supported by two other observations. First, its position is consistent with the 
theoretical expectations — moving from 1/4 at strong coupling to 1/8 at weak coupling. Second, 
the jump at the transition decreases as b increases, consistent with the theoretical expectation that 
the transition is driven by lattice artefacts. 
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4. Conclusions and Outlook 

In summary, our evidence to date is consistent with the conjecture in the right panel of Fig. |I|, 
including the presence of a window around jq. of finite width in which the symmetry is un- 
broken. Accepting this result, the adjoint Nf = I theory provides the first successful example of 
reduction to a single-site. Furthermore, it appears that the funnel region is broad enough to allow 
one to work with quarks with masses m » Aqcd (and possibly with m I /a) so that the long- 
distance physics is that of the pure gauge theory. In this case we have, in effect, a version of the 
original EK model in which the effect of the adjoint quarks is to maintain the repulsion among 
eigenvalues of the Polyakov loop, presumably by inducing a tower of double-trace operators simi- 



lar to that suggested in Ref. [ ]14[ ] as a cure for the EK model. If, on the other hand, we work near 
Kc we have a single site model that lies "within 1 /N" of physical two-flavor QCD. 

Our analysis can be extended in several ways. Our interpretation of the critical line can be 
checked by calculating the pion mass. We can also calculate other meson masses and the string 
tension, and the distribution of the eigenvalues of Dirac operators for various valence fermions. 
Work in these directions is underway using the configurations in hand. We also plan on speeding 
up the algorithm so as to move to larger values of N. 

Finally, we note that another recent lattice study (using overlap as opposed to our Wilson 



fermions) finds preliminary evidence for a failure of reduction QISQ. This discrepancy with our 
findings will clearly need to be resolved. 
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